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Abstract: This work presents preliminary results for the time-dependent cosmic-ray propagation in the 
Galaxy by a fully 3-dimensional Monte Carlo simulation. The distribution of cosmic-rays (both protons 
and helium nuclei) in the Galaxy is studied on various spatial scales for both constant and variable cosmic- 
ray sources. The continuous diffuse gamma-ray emission produced by cosmic-rays during the propagation 
is evaluated. 



Introduction 

One of the important goals of current cosmic- 
ray studies is to understand the origin as well as 
the propagation of cosmic-rays, in particular the 
sources of cosmic -ray particles and also the distri- 
bution of sources in the Galaxy. One widely used 
method is to examine the matter path length, or 
age distribution of cosmic -rays [2, 8]. Originally, 
the most popular model to describe the cosmic- 
ray propagation in the Galaxy was the leaky-box 
model (LBM) [4] that unrealistically assumes a 
uniform source distribution and the diffusion and 
boundary effect are taken into account by means 
of simple loss terms. Diffusion models [1] are 
more realistic but their analytic forms are limited 
in the description of the many parameters required 
and their spatial and energy dependence that are 
needed to fully describe cosmic-ray propagation in 
the Galaxy. A numerical solution of the diffusion 
equation by a finite-difference scheme avoids this 
problem, but is only accurate if the grid is much 
finer than all physical scales [9]. A Monte Carlo 
approach should be ideal for the study of cosmic- 
ray propagation in the vicinity of their sources. 

This work presents preliminary results for the 
time-dependent cosmic-ray propagation in the 
Galaxy by a fully 3-dimensional Monte Carlo sim- 
ulation. The distribution of cosmic -rays in the 
Galaxy is studied on various spatial scales for both 
constant and variable cosmic-ray sources. The 



continuous diffuse gamma-ray emission during the 
cosmic-ray propagation is evaluated. 

Model and Simulation Technique 

In the simulation, the cosmic-ray diffusion is de- 
scribed by a Monte Carlo process in which the par- 
ticle position at time t n +\ is related to its position 
at t n by 

Xi(t n+ i) = Xi(t n ) + Axi = Xi(t n ) + cosaiAr (1) 

where the particle 3-coordinate, Axi the dis- 
placement for each coordinate, cos on the direction 
cosine for each direction determined by a random 
number generators and Ar the mean displacement 
that is described as 



Ar = y/6D(E)At 



(2) 



with D(E) as the energy-dependent diffusion co- 
efficient and the timestep At. The choice of At 
directly determines the spatial resolution. 

The simulation methodology for the particle mo- 
tion through physical space and energy space is 
as follows. After reading the physical configura- 
tions such as Galactic dimension, ISM distribution, 
Galactic magnetic field, total propagation time, de- 
tection shell dimension, source dimension and in- 
jection profile (energy range, spectral profile, etc.), 
a particle, after been injected by a source, is al- 
lowed to undergo a fully 3-dimensional random 
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walk. At each step, the particle experiences en- 
ergy loss and interactions. The particle contin- 
ues to propagate until either it reaches the bound- 
ary, or the propagation time scale is greater than 
10 8 years, or its energy is below the minimum sim- 
ulation energy (i.e., -Emi n = m p + 0.3 GeV ~ 
1.24 GeV). After reaching the end of the simu- 
lated propagation of a particle, the next particle is 
simulated. This technique allows the parallel com- 
putation without communication. 

Result Examples from Simulation 

The transport equation of cosmic rays is generally 
written as 



dN 

~dt 
_d_ 

dE 



Q(E,r,t) + W ■ (DWNi) - Pi Ni 
[b i (E)N l (E)] — V • (VNi(E)) 



k>i 



where Ni(E, r, t)dE is the particle number den- 
sity of type i at position r with energy in (E, E + 
dE). Other terms, e.g. for stochastic reaccelera- 
tion, can be added as desired. For a simple test of 
the simulation code, we simplify Eq. (3) by keep- 
ing only the source term and the diffusion term: 



dN 
~dt 



Q{E, r, t) + V • (DVJV). 



(4) 



Without losing generality, the solution to Eq. (4) 
can be written as 



N(E,r,t) 



dV' dE' dt'Q-G (5) 



by the Green's function such that 



G(E,r,t) 



1 



8(nDt) 3 / 2 



exp 



7(40*)] (6) 



where G(r,t) gives the probability of finding a 
particle at position r and at time t that was injected 
at r' = and t' = 0. 

Fig. 1 shows the simulation results (dashed lines) 
and the analytic solutions (solid lines) calculated 
by Eqs. (5) and (6) for the particle number in dif- 
ferent shells, following the instantaneous injection 
of 10 5 particles at the origin. 



o 

£ 

o 
o 




10 10 

Elapsed Time (Yr) 

Figure 1: Analytic solutions (solid lines) and 
simulation results (dashed lines) for the diffu- 
sion equation (4). We count particles in shells 
of width 1 pc, that are located at |r| = 
10, 50, 100, 500 and 1000 pc. A total of 10 5 par- 
ticles are simulated. 

We also consider a case with continuous energy 
loss by synchrotron radiation and inverse Comp- 
ton scattering, and an energy-dependent diffu- 
sion coefficient D(E) = DqE s for the case 
of Kolmogorov-type spectrum of turbulence [5]. 
Thus the transport equation can be reduced to 

dN d 

~dt~dE ^ E2N) - d oE S V 2 N^Q (7) 
with the source term 

Q = Q Q E- a O(T-t )e(t + T-T). (8) 
The analytic solution to Eq. (7) is obtained as [7] 



N = Q Q E~ a / dr 

J-l/b E 

8(t - t )e(t +T-t) cxp(-r 2 /4A) 
(47tA) 3 / 2 (1 + 6 £t) 2 ~ q 

where the cooling time r and the length scale A are 

f E " dE' f E " D(E')dE' 

Fig. 2 shows the simulation results (dashed lines) 
and the analytic solutions (solid lines) to Eq. (7) 
with an injection of total 10 5 particles with an ini- 
tial energy E — 200 GeV by a source at t — and 
at the origin. The figure is arranged as Fig. 1. 



(9) 
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Figure 2: Analytic solutions (solid lines) and simu- 
lation results (dashed lines) for the diffusion equa- 
tion (7) for 10 5 particles with a same energy in- 
jected by the source. The figure is arranged as 
Fig. 1 but with shell width of 10% of \r \ for each 
shell for better statistics. 

Finally we consider a power-law source spectrum 
Q = QoE~ a for cosmic-ray ions. To simu- 
late with the same statistical accuracy at all ener- 
gies, each simulated particle carries a weight factor 
which is calculated for each energy bin: 

w (E) = Q E- a ^SE (11) 

and the normalization factor Q is obtained by 

J QdE = J Q E~ a dE = N (12) 

with Nq the total particle number in the simulation 
and A^Bin the total energy bin number. Fig. 3 shows 
the simulation results (dashed lines) and analytic 
solutions (solid lines) of the particle distribution 
in shells 1 and 2 at different elapsed time: in the 
ranges of 10 2 — 10 4 years and 10 4 — 10 5 years, 
respectively, after been instantaneously injected by 
a source at the origin and with an inelastic energy 
loss . Note a total of 10 7 particles with energies 
from 10 to 10 5 GeV are simulated in this case. The 
inelastic energy loss in this energy range is approx- 
imated by E ~ 5 • lO- 16 ^ 1 - 07 [6] to facilitate the 
analytical solution and replace the synchrotron en- 
ergy loss rate in Eqs. (9) and (10). Also note the 
time is relatively short when compared with the 
particle propagation time (10 8 years), therefore we 
describe the particle distribution in the surround- 
ings of sources. 



Diffuse 7-Ray Production 

By using a 7-ray production matrix published ear- 
lier [6], the diffuse 7-ray spectrum produced by the 
cosmic-rays can be determined. Fig. 4 shows the 
continuous 7-ray spectra produced by cosmic-rays 
with distributions at different elapsed timescales 
shown in Fig. 3. Figs. 3 and 4 present an exam- 
ple of spatial distribution and time evolution of 
cosmic-rays and 7-rays in the source vicinity. The 
cosmic-ray population in the Galaxy is thought 
to be approximately stationary, even though occa- 
sional spikes in the cosmic -ray flux must be ex- 
pected [3]. This study of cosmic -ray propagation 
may provide a more detailed view of spectral vari- 
ations in the vicinity of cosmic -ray sources. 

Conclusions 

Cosmic-ray propagation in Galaxy and 7-ray pro- 
duction during the cosmic-ray propagation are 
studied by a fully 3-dimensional Monte Carlo sim- 
ulation approach. Our results are in good agree- 
ment with analytic solutions to the diffusion equa- 
tion. The spatial resolution in this simulation 
varies between about 0.5 pc and about 10 pc, de- 
pending on particle energy as well as the elapsed 
time. The calculations may prove useful in de- 
termining the distribution and spectra of cosmic 
rays and their time evolution in the vicinity of their 
sources. 
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Figure 3: Particle distribution calculated by simulation (dashed lines) and analytic solutions (solid lines) for 
a power-law injection by a source at the origin and at t = 0. Figure shows examples for 2 shells as defined 
in Fig. 1 but with shell widths 20% of r for different elapsed time: in range of 10 2 -10 4 years (shell 1) and 
10 4 -10 5 years (shell 2). 
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Figure 4: Continuous 7-ray energy spectra produced for a power-law cosmic -ray spectrum injected by a 
source at the origin and at t = 0. Figure shows results for shells 1 & 2 as defined in Fig. 3. 
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